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Abstract 



We construct a bound state of three 1/3-quantized Josephson coupled vortices in 
three-component superconductors with intrinsic Josephson couplings, which may be rel- 
evant with regard to iron-based superconductors. We find a Y-shaped junction of three 
domain walls connecting the three vortices, resembling the baryonic bound state of three 
quarks in QCD. The appearance of the Y-junction (but not a A-junction) implies that 
in both cases of superconductors and QCD, the bound state is described by a genuine 
three-body interaction (but not by the sum of two-body interactions). We also discuss a 
confinement /deconfinement phase transition. 
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1 Introduction 



Color confinement in quantum chromo dynamics (QCD) is one of the most challenging problems in 
modern physics. Quarks having fractional electric charges ±(l/3)e or ±(2/3)e should be confined 
by color electric fluxes to form hadrons (mesons or baryons) with integer electric charges. As a 
result, only mesons and baryons, made of two and three quarks, respectively, can be observed. A 
color electric flux tube stretched between a quark and an anti-quark provides constant attractive 
force or a potential that is linearly dependent on the distance between the quarks. A recent 
lattice QCD simulation has confirmed this picture [lj. On the other hand, one can imagine two 
possible configurations of color fluxes in a baryon: a Y-shaped junction or a A-shaped junction 
connecting three quarks. The A-junction implies that the three-quark interaction is described by 
the sum of two-body interactions while the Y-junction implies a genuine three-body interaction. 
Although it has been a long standing issue as to which of these possible configurations is the 
actual one, a further study of lattice QCD simulations has clearly demonstrated the Y-junction 

A bound state of Josephson coupled vortices confined by domain walls exists in multicom- 
ponent superconductors, and it resembles the bound states of quarks confined by fluxes. Half- 
quantized vortices stably exist in two-band superconductors when the interband Josephson cou- 
pling is negligible [3J. When the interband Josephson coupling is taken into account, a half- 
quantized vortex is attached by a domain wall [U [3] which extends to the edge of the sample. 
Since the domain wall pulls the vortex in order to reduce the energy, the half-quantized vortex 
is unstable. On the other hand, two half-quantized vortices winding around two different gap 
functions are connected by a domain wall (sine-Gordon kink) [U |5] when the interband Josephson 
coupling is considered. The domain wall provides an attractive potential linearly depending on 
the distance between the two vortices [6j; this bound state resembles a meson in QCD. In fact, 
confinement /deconfinement phase transition occurs at finite temperature, similar to the case in 
QCD, in which fluctuations of the domain wall contribute to entropy [6J. That is, vortices with 
fractional quanta are confined to become vortices with integer quanta below a certain tempera- 
ture, while they are deconfined above that temperature. However, the question is: Is there any 
model or actual material that indicates a baryonic bound state of vortices? In the absence of 
intrinsic Josephson terms, vortex bound states have been discussed in multi-component super- 
conductors [7]. How are they connected by domain walls, once the Josephson terms are turned 
on? 

In this paper, we show that baryonic bound states indeed exist in three-component supercon- 
ductors with intrinsic Josephson terms. We show that three 1/3-quantized vortices are connected 
by a Y-shaped junction of domain walls, resembling a baryon in QCD. We also discuss a confine- 
ment /deconfinement phase transition. Above the certain critical temperatre, the three fractional 
vortices are deconfined from a vortex baryon. 

Some two-band superconductors exhibit type 1.5 superconductivity, i.e., repulsion at a short 
distance and attraction at a large distance between two integer vortices j8]. Such a structure 
leads to a cluster of vortices, which has been experimentally confirmed in MgB2 [HI [10]. Recently, 
three-component superconductors have attracted considerable attention because of the discovery 
of iron-based superconductors [UJ. Therefore, we propose the multi-band superconductors to test 
baryonic bound states of vortices with 1/3 quanta which exhibit a confinement /deconfinement 
transition. Our solution will also suggest the possibility of experimentally determining intrinsic 
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Josephson couplings by determining the shape (angles and lengths) of a Y-junction of three 
vortices. 



2 Three-component superconductors 



Multicomponent (n-component) superconductors can be generically described by the Ginzburg- 
Landau free energy, 
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with the intrinsic Josephson terms 
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Here, 7^ are constants and is decomposed into amplitude and phase as = \^i\e \ The 
phases of and ^ preferably coincide for 7^ > while they tend to have tt phase difference 
for jij < 0. All phases are the same in the ground state when 7^ > for all % and j, while the 
system is frustrated when 7^ < for all i and j. The Hamiltonian (CQ) is invariant under gauge 
symmetry, 



he 
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In the limit of 7^ = 0, the Hamiltonian (JTJ) enjoys t/(l) n symmetry, of which the overall phase 
rotation exhibits the gauge symmetry ()3]) while others exhibit global symmetry. In this case, 
there appear n — 1 Nambu-Goldstone modes. They are gapped for non-zero 7^ (the Legget 
modes). 

Hereafter, we restrict ourselves to three-component (n = 3) superconductors in two space 
dimensions [T2J[32]- We consider the positive Josephson couplings, 7^ > 0, while the negagive 
Josephson couplings, 7^ < 0, give frustrated systems [T3]. There exist three types of vortices, 
labeled as (1, 0, 0), (0, 1, 0), and (0, 0, 1), winding around the first, second, and third component by 
27r, respectively. The energy of each vortex is logarithmically divergent when 7^ = and linearly 
divergent when 7^ 7^ 0, if the system size is infinite. Let us consider that the (1,0,0), (0, 1,0), 
and (0, 0, 1) vortices are placed at the edges (Pi, P2, and P3, respectively) of a triangle, as shown 
in Fig. [TJ In the large circle, the total configuration is the integer vortex (1, 1, 1), which implies 
that the total energy is finite. In other words, the integer vortex has an internal structure made 
of three vortices, all of which are 1/3 quantized, as shown below. 

Instead of the U(l) 3 generators (1, 0, 0), (0, 1, 0), and (0,0,1), let us prepare four linearly 
dependent generators: the gauge rotation (1, 1, 1) and three gauge-invariant rotations (0, —1, 1), 
(1, 0, — 1), and (—1, 1, 0). Among these, only the gauge rotation is accompanied by gauge trans- 
formation ()3]), while the others are all global phase rotations. In these new generators, the 



2 




Figure 1: (Color online) The (1, 0, 0), (0, 1, 0), and (0, 0, 1) vortices are placed at Pi, P2 and P3, 
respectively. 6j (i = 1, 2, 3) corresponds to 1/3 circles at the boundary, and corresponds to the 
radial paths from the origin O to the circle at the boundary The (1, 0, 0), (0, 1, 0), and (0, 0, 1) 
vortices are encircled by 61 — r 3 + r 2 , b 2 — 77 + r 3 , and 63 — r 2 + 77, respectively. 

winding of the (1, 0, 0), (0, 1, 0), and (0, 0, 1) vortices can be decomposed into 

Pi : (1, 0, 0) = ^(1, 1, 1) + 0(0, -1, 1) + ^(1, 0, -1) - ^(-1, 1, 0), 
P 2 : (0, 1, 0) = 1, 1) - i(0, -1, 1) + 0(1, 0, -1) + 1, 0), 

P 3 : (0, 0, 1) = 1(1, 1, 1) + ~(0, -1, 1) - 0, -1) + 0(-l, 1, 0). (4) 

We see that all the paths b{ (i = 1,2,3) in Fig. [TJ correspond to 27r/3 rotation of the gauge 
generator (1, 1, 1) with Eq. ([3]) and consequently that these vortices are all 1/3 quantized; their 
magnetic flux is $o/3 with the unit flux quanta $0 = hc/2e [H]. The (1,0,0), (0,1,0), and 
(0, 0, 1) vortices are encircled by b\ — r 3 + r 2 , 62 — 77 + r 3 , and 63 — r 2 + 77, respectively (Fig. [TJ). 
Therefore, we can identify the paths ±77, ±r 2 and ±r 3 corresponding to ±27i~/3 of the global 
phase rotations (0, —1, 1), (1, 0, —1), and (—1, 1, 0), respectively. The global phase rotation along 
the radial path r,i can be written up to constant phases as 

77 : ^ = 1^1, ^ 2 = e- (2ni/3Wr) \y 2 \, ^3 = e {2 ™ /3)/(r) |tf 3 |, 
r 2 : *j = e^/^/Wi^^ ^ = ^ 3 = e -(2«/3)/(r) 

r 3 : = e - (2 " /3)/(r) |*i|, * 2 = e (2wi/3)/(r) |^ 2 |, * 3 =|* 3 |, (5) 

where a function /(r) has the boundary conditions f(r = 0) = 1 and f(r — > 00) = 0. 

The integer vortex configuration (1,1,1) is of the Abrikosv type which has finite energy (due 
to the fact that the associated U(l) symmetry is fully local), whereas the fractional vortices cor- 
respond to global ungauged symmetries and hence they have a logarithmically divergent energy, 
even in the absence of the Josephson terms, 7^ = 0. This is why the fractional vortices are 
confined whereas the integer vortices are acceptable finite-energy solutions. 

In the presence of the Josephson terms, 7^ 7^ 0, we expect there to be a sine-Gordon kink 
in each path r, that connects two vortices. However, the question that arises is: how does it 
connect two vortices? Does it connect along the segment PjPfc? 
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3 Baryonic bound state 



We concentrate on the case of 7^ > for all i and j, which may be the case of iron-based 
superconductors. In this case, all phases are the same in the ground state, which is unique 
with respect to gauge transformation ([3]). The phases at r — > 00 are the same, which we 
consider to be zero according to the gauge symmetry; hence, Eq. <^ holds, including for the 
constant phases. Along each radial path, the Josephson term in Eq. (J2]) can be written as 
7ij|\l/j| cos((47r/3)/(r)), from Eq. ([5]). We thus find that it takes a non-zero value oTy'l^Ml^il 
at the center (r = 0). Therefore, the linear connection of vortices along the segments PjPfc due 
to the sine-Gordon kinks would increase energy because a domain (membrane) with finite energy 
appears inside the triangle P1P2P3. The sine-Gordon kinks should bend to form the Y-junction. 

In fact, the numerical solution indicates the Y-junction, as seen in Fig. [2j For simplicity, we 
have taken 7^ = 7 (> 0) but the general case is straightforward. In Fig. [2]-(c), we find that the 
magnetic field is localized at the center of each vortex. Fig. [3]-(a) shows the phases of ^2, 
and \l/3, indicating a phase winding at Pj. We also obtain the numerical solution for the function 
f(r) in Eq. ([5]), (Fig. [3]-(b)). In this numerical simulation, we have fixed the positions of the 
three vortices. The wall tension leads to a linear potential (confining force), and these vortices 
collapse to form a single integer vortex. 



4 Confinement /deconfinement phase transition 

Here we discuss that the Y-junction can be stable at finite temperature and exhibits a phase 
transition, as proposed by Goryo et.al. [6 J for a mesonic bound state of vortices in two-gap 
superconductors. To this end, we construct the effective theory for the phases of the gap functions 
with keeping the amplitudes constants, ^/«( x ) = v i ex P(^«( x ))) given by 



eff. 



(6) 



up to a constant. Let us define the phase diferences as gauge invariant dynamical variables as 

^ = 03 - 03, fo = 9 3 -9 1 , 03 = 0! -0 2 , (0, + 2 + 03 = 0). (7) 

With taking a gauge fixing as 

0i + 2 + 3 = 0, (8) 

the energy (jBJ) is further reduced to 

~h 2 v? 



"cS. 



E 



12m 



(V&r + 77f(l-cos. 



(9) 



where we have set A = and defined 



Vl = ^3723, V2 = U3«l731, = VlV 2 -)l2- 



(10) 
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(b) (c) 

Figure 2: (Color online) Baryonic bound state of vortices. Plots of (a) the total energy density, 
(b) the energy density of the Josephson couplings, and (c) the magnetic field. For simplicity, we 
take jij = 7 > 0, but the general case is straightforward. The relaxation method has been used 
with parameters h = c = 2m = 2e = v = A/2 = 1 and 7 = 0.02. 




Figure 3: (Color online) (a) The arrows and contour lines indicate, respectively, the phases and 
amplitudes of (left) , ^2 (middle) and ^3 (right), (b) The plot of the function f(r). 



5 




(a) (b) 

Figure 4: (a) The two vortices P 2 (0,1,0) and P 3 (0,0,1) together are placed at the same 
position very far from the vortex Pi (1, 0, 0). They are connected by a sine-Gordon domain wall, 
(b) The most symmetric configuration at a finite temperature. 



In order to calculate the tension of the domain wall attached to the vortex (1,0,0), let us 
place the two vortices P 2 (0, 1, 0) and P3 (0, 0, 1) together at the same position which is very far 
from the vortex P 1 (1, 0, 0) as in Fig. I3ja). In this situation, we can set 9 2 = 63 so that we have 

01 = 0, 2 = -03 = 0- (11) 

Then, the effective model ([HD reduces to the sine-Gordon model 

F«=K«(V0) 2 + ^(l-cos0), 

K (l)- V(vl + Vl) r (l)_ 2( 2 2) (12) 

12m ' W2 l3J K ' 

This can be rewritten as the Bogomol'nyi form 

F&] = [^ 1 )(V0) 2 + r( 1 )sin 2 (0/2)] 

= (VKWV(f) ± VfW sin(0/2)) 2 =f 2v / K( 1 )r( 1 )V0sin(0/2) 

> e (13) 

with the topological charge density in the second term in the second line, 

e = ^2v / KWT« V0 sin(0/2) = ±AVKWTWV cos(0/2). (14) 

The most stable configurations with the minimum energy can be achieved by satisfying the 
Bogomol'nyi equation, obtained by (....) 2 = in the second line of Eq. f[T3"|) . i.e. 

v^WV0± v / f«sin(0/2) = 0. (15) 
One (anti-)kink solution can be obtained as 



4 arctan exp 



+ 1 / r(1) f \ 



(16) 
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where x is the coordinate perpendicular to the kink and xq denotes the position of the kink. The 
tension of the one (anti-)kink is 



dx e 



WKWTW |[cos(0/2)]*I+~| = 



h 2 (v 2 2 + vl){r,l + rj) _ vi(vj + vj)(v2^i + vslsi) {m 
12m V 6m ' 1 ' 

In the same way, the other two domain walls attached to the (0, 1,0) and (0,0, 1) have the 
tensions 

= 8VK(W*>, = h2 ( v l + v D r (2) = 2(2+ 2) 

12m ' W3 llh 

T( 3 ) = sVPlrW, = ^M±j|) } r ( 3 ) = 2(r /l 2 + r /2 2 ), (18) 

respectively. 

We are now ready to discuss the confinement /deconfmement phase transition. For simplicity, 
we consider the most symmetric case with 

v 1 = v 2 = v 3 = v, 7i 2 = 723 = 731 = 7- (19) 
The tension of each domain wall becomes 



T dw = 8^ 2 ./^. (20) 
V om 

In this case, the molecule is Z 3 symmetric as in Fig. H^b). For each domain wall with the length 
L, the total energy and entropy can be evaluated as [SI EH] 

E = T dw L, 5*dw = t L, (21) 

respectively, with a short lengh cut-off £ which is the largest among the coherence length and 
the penetration depth. Consequently, the free energy of each domain wall at the temperature T 
is given by 

ksT In 2-7T 



-fdw — Edw — TSdw — ( Tdw 1 J L — AL, (22) 

A = T dw -^^ . (23) 

When the coefficient A is positive, the interger vortex is stable, i.e., in the confinement phase. 
On the other hand, when the coefficient A is negative, the integer vortex tends to be split into 
a set of the three fractional vortices in order to reduce the free energy of the domain walls, that 
is, the deconfinement occurs. Therefore, the critical temperature for the deconfmement is found 
to be 



Tcrit = ^ = l ^ hv / 27 . (24) 
cn /csln27r A;£ln27rV3m 

This expression is the same with Goryo et. al [6] . In the most symmetric case with Eq. ( TT9|) which 
we are considering, the confinement mechanism is essentially the same with the case of two-gap 
superconductors. 
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5 Summary and Discussion 



In summary, we have constructed baryonic states of three 1/3-quantized vortices and have found 
that these vortices are connected by the Y-junction of domain walls, resembling a baryonic bound 
state of three quarks in QCD. In both cases of superconductors and QCD, the appearance of the 
Y-junction and not of a A-junction implies that the bound state is described by a genuine three- 
body interaction and not by the sum of two-body interactions. The confinement /deconfmement 
transition of vortices has been studied. This common feature between superconductors and QCD 
will shed a new light on the color confinement problem of QCD. 

The similarities between superconductors and QCD should be further clarified. As a toy model 
of QCD, Shifman and Unsal [19] considered an SU(2) gauge theory in three space dimensions 
with one direction compactified as R 2 x S 1 . The theory becomes a U(l) gauge theory in two 
dimensional space R 2 in the limit of a small radius of S 1 . It was shown by Polyakov [20J that the 
confinement occurs in a U(l) gauge theory in two dimensional space; electrically charged particles 
(quarks) are confined by an electric flux. By taking a duality, quarks are mapped into vortices, 
while electric fluxes are mapped into sine-Gordon domain walls, so that a meson made of two 
quarks is mapped to a mesonic bound state of two vortices. There, the quark confinement can 
be understood as the vortex confinement which is described by two-gap superconductors. The 
confinement is nontrivial in the former, while it can be easily shown in the latter. We therefore 
expect that a discussion along the same line shows that a baryon made of three quarks in SU(3) 
QCD is mapped to a bayonic bound state of three vortices found in this paper. We conjecture 
that the existence of a Y-junction in a baryon in QCD can be shown by using a duality map to 
three-gap superconductors. In multi-gap superconductors with more than three gaps, the bound 
states of more vortices should exist which may correspond to tetraquarks, pentaquarks, etc. in 
QCD. 

Although we have chosen the same Josephson couplings 7^ = 7 as an example, an extension to 
the general case is straightforward. It may be useful to determine intrinsic Josephson couplings of 
multicomponent superconductors such as iron-based superconductors by determining the shape 
of the Y-junction. Several interesting studies have been conducted on for vortex mesons in two- 
component or p-wave superconductors, for instance, the studies on a lattice of vortex mesons and 
twistons [21] and those on vortex clusters [9j HD] • We hope that our work will stimulate further 
theoretical and experimental studies of multicomponent superconductors, particularly iron-based 
ones. 
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